Detection of cosmic superstrings by geodesic test particle motion 
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(p,q)-strings are bound states of p F-strings and q D-strings and are predicted to form at the 
end of brane inflation. As such these cosmic superstrings should be detectable in the universe. In 
this paper we argue that they can be detected by the way that massive and massless test particles 
move in the space-time of these cosmic superstrings, in particular we study solutions to the geodesic 
equation in the space-time of field theoretical (p,q)-strings. The geodesies can be classified according 
to the test particle's energy, angular momentum and momentum in the direction of the string axis. 
We discuss how the change of the magnetic fluxes, the ratio between the symmetry breaking scale 
and the Planck mass, the Higgs to gauge boson mass ratios and the binding between the F- and 
D-strings, respectively, influence the motion of the test particles. While massless test particles can 
only move on escape orbits, a new feature as compared to the infinitely thin string limit is the 
existence of bound orbits for massive test particles. In particular, we observe that - in contrast to 
the space-time of a single Abelian-Higgs string - bound orbits for massive test particles in (p,q)- 
string space-times are possible if the Higgs boson mass is larger than the gauge boson mass. We 
also compute the effect of the binding between the p- and the q-string on observables such as the 
light deflection and the perihelion shift. While light deflection can also be caused by other matter 
distributions, the possibility of a negative perihelion shift seems to be a feature of finite width cosmic 
strings that could lead to the unmistakable identification of such objects. In Melvin space-times, 
which are asymptotically non-conical, massive test particles have to move on bound orbits, while 
massless test particles can only escape to infinity if their angular momentum vanishes. 
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I. INTRODUCTION 



Cosmic strings are topological defects that are predicted to have formed via the Kibble mechanism [l| during one 
of the phase transitions in the early universe and in the field theoretical description 0] can be considered to be an 
example of a topological soliton. Due to the fact that these objects can be extremely heavy (up to 10^^ kg/m) they 
were believed to be a possible source of the density perturbations that led to structure formation and the anisotropics 
in the cosmic microwave background (CMB) [3]. However, the detailed measurement of the CMB power spectrum 
as obtained by COBE, BOOMERanG and WMAP showed that cosmic strings cannot be the main source for these 
anisotropics. 

In recent years cosmic strings gained renewed interest due to the possible connection to the fundamental entities 
of String Theory Q . Brane inflation is a popular inflationary model that can be embedded into String Theory and 
predicts the formation of cosmic string networks at the end of inflation 0]. E.g. in the framework of type IIB 
String Theory the inflaton field corresponds to the distance between two Dirichlet branes with 3 spatial dimensions 
(D3-branes) and inflation ends when these two branes collide and annihilate. The production of strings (and lower 
dimensional branes) then results from the collision of these two branes. Each of the original D3-branes has a U(l) 
gauge symmetry that gets broken when the branes annihilate. If the gauge combination is Higgsed, magnetic flux 
tubes of this gauge field carrying Ramond-Ramond (R-R) charge are D-branes with one spatial dimension, so-called 
D-strings. When the gauge combination is confined the field is condensated into electric flux tubes carrying Neveu 
Schwarz-Neveu Schwarz (NS-NS) charges^and these objects are fundamental strings (F-strings) [6]. D-strings and F- 
strings are so-called cosmic superstrings [4|] which seem to be a generic prediction of supersymmetric hybrid inflation 
and grand unified based inflationary models Q. D- and F-strings, however, have different properties than the 
usual (solitonic) cosmic strings. The probability of intercommutation of solitonic strings is equal to one but less than 
one in the case of cosmic superstring. Therefore, solitonic strings do not merge, while cosmic superstrings tend to 
form bound states. When p F-strings and q D-strings interact, they can merge and form bound states, so-called 
(p,q)-strings whose properties have been investigated [l^l- Even though the origin of (p,q]-strings is type IIB 
string theory, their properties can be investigated in the framework of fleld theoretical models [Ill4l4| . The influence 
of gravity on field theoretical (p,q)-strings has been studied in (isj . 

Since there are good reasons for cosmic superstrings to be a consequence of String Theory it is very exciting to 
search for observational consequences of their existence. There has been considerable effort in numerically modeling 
cosmic string networks to obtain CMB power and polarization spectra . Comparison with observations has shown 
that cosmic strings might well contribute considerably to the energy density of the universe. There is another way to 
detect cosmic strings in the universe, namely through the motion of test bodies in such strin g sp ace-times. The test 
particle motion in different space-times containing cosmic strings has been investigated in |17h21| , while the complete 
set of orbits of test particles in the space-time of black hole pierced by an infinitely thin cosmic string has been given 
for a Schwarzschild black hole in [23] and for a Kerr black hole in [23]. 

In this paper we follow the latter approach and use the field theoretical model discussed in [l^ to describe (p,q)- 
strings by two coupled Abelian-Higgs models in curved space-time. For vanishing coupling between the two sectors, 
the model corresponds to the Abelian-Higgs model coupled minimally to gravity. This model has solutions describing 
strings with finite core width that have been investigated in 0, [2^ . Geodesies in this space-time have been studied 
recently 26] and can only be given numerically. Here we would like to extend this investigation to the field theoretical 
description of cosmic superstrings. 

Our paper is organized as follows: in Section II, we discuss the field theoretical model that possesses (p,q)-string 
solutions and we also work out the geodesic equation. In Section HI we discuss our numerical results, in particular we 
give examples of orbits and demonstrate how the ratio between the symmetry breaking scale and the Planck mass, 
the ratios between Higgs and gauge boson masses, the magnetic fiuxes and the binding between the F- and D-string 
infiuence our results. We conclude in Section IV. 



II. THE MODEL 
A. The space-time of a (p,q)-string 

The field theoretical model to describe gravitating (p,q)-strings reads [iB] 



I d''xy^[-^R + C^,] , (1) 



IGttG 
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where R is the Ricci scalar and G is Newton's constant. The matter Lagrangian is given by 

Crn = D,cf,{D^^c^r ~~ iF^.i^'"- + D^^D'-O* - \h^,H'^-' ~ u{<j>, (2) 

with the covariant derivatives Df^(j> — - ieiA^cf), D^£^ — V^^ - ie-zBf^S, of the two complex scalar fields (Higgs fields) 
(f) and ^ and the field strength tensors Ff^i, = V^Ai. - Vi.A^ = d^A^ - di,A^, H^^y = V ^Bi, ~ Vi,B^ = d^Bi, - di,B^ of 
two U(l) gauge potential A^, B^, with coupling constants ei and 62- denotes the gravitational covariant derivative. 
Finally, the potential V{(j),(,) reads: 

"(0,0 = ^i^r iif + ^(^r - riif - A3(00* - lim* vD , (3) 

where Ai and A2 are the self-couplings of the two scalar fields, while A3 > is the coupling between the two sectors, 
ryi and 772 are the vacuum expectation values of the scalar fields. 

In order for both U(l) symmetries to spontaneously break which then leads to the formation of strings we have to 
require that the (absolute) minimum of the potential ([3]) is at non- vanishing values of |(/)| and |^|. This leads to the 
requirement 

A1A2 > A\l . (4) 
The most general static cylindrically symmetric line element invariant under boosts along the z-direction is 

ds^ = N{pfdf - _ Lipfdi/ - N{pfdz^ . (5) 
For the matter and gauge fields, we apply the Ansatz [2] 

0(p,^) = r;i;i(p)e™^ , A^dx^ = -{n-P{p))d^ (6) 

ei 

^(p,V') = ^2/(p)e™'^ , B.dx^^ ^-{m-R{p))d^ , (7) 

62 

where n and m are integers indexing the vorticity of the two Higgs fields around the z-axis and correspond to the 
degree of the map from — >■ S^, where the homotopy group is 7ri(S'^) = Z. In our field theoretical model of 
(p,q)-strings the p corresponds to the winding n and the q to the winding m. 
We can then do the following rescaling 

p^-^ , L^— (8) 
eiTyi eiTji 

such that the total Lagrangian only depends on the following dimensionless coupling constants 

Q 2 62 772 o fn\ 

ei T]i ef 

where i = 1,2,3. 7 is proportional to the ratio between the Planck mass A/pi = G^^^^ and the symmetry breaking 
scale iji. Moreover, ^/pi is proportional to the ratio between the Higgs mass Mh,i = V^Vi and the corresponding 
gauge boson mass Afyv,! — V^eirji, while \/^/g is proportional to the ratio between the Higgs mass Mh,2 = \fy^f]2 
and the corresponding gauge boson mass M-^^i — A/2e2'72- Each of the strings possesses a scalar core with width 
PH.i ~ -^^H i ^^"^ ^ gauge field core with width pw.i ~ ^^wY i = 1, 2. Note that with the rescaling ([8|) the width of 
the gauge field cores is pw.i ^ l/V^, Pw,2 ~ \l{gq\f2) while the widths of the scalar cores is given by pH.i = l/^/Pl, 
i = l,2. 

The variation of the action ([ij with respect to the matter fields leads to the following equations [l^ 

(10) 
(11) 
(12) 

(13) 





{N^Lh'Y 


P^h 1 du 




N^L 
(N^LfJ 


+ 2dh 
R^f ^ Idu 






^ 2df 


L 


rn 


= 2h'^P , 


L 


(77 


f 
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where the prime denotes the derivative with respect to p and the potential u reads 



/3; 



The variation of ([1]) with respect to the metric leads to the Einstein equations 
where T is the trace of the energy momentum tensor. Using our Ansatz these read [l£ 



'2/ .2 



{LNNJ 

{NHjy_ 



{P'f [R'f 



2L2 2g2L2 

2/i2p2 2i?2/2 (p')2 



{R' 



L2 



L2 



2L2 2g2L2 



In addition there is a constraint equation that is not independent. This reads 



' NL 



{N 



l\2 



iV2 



7 



7\2 



'^2 



/l2p2 



2i2 252L2 £2 



2 f2 



L2 



(14) 
(15) 

(16) 
(17) 

(18) 



The set of differential equations can be solved only numerically subject to an appropriate set of boundary conditions. 
The requirement of regularity at p = leads to the following conditions 



/i(0) = 0, /(O) = 0, P(0) = n, P(0) = m 

for the matter fields and 

N{Q) = 1, A^'(O) = 0, i(0) = 0, L'(0) = 1 
for the metric fields, while the requirement of finiteness of the energy per unit length leads to 

h{oo) = 1, /(oo) ^ q , P(oo) = , P(oo) = . 
The inertial energy per unit length of the (p,q)-string is given by 



2tt 



NL {h'Y + {f'Y + 



(P' 



{R'Y h^P^ P2j2 



2L2 2g2L2 12 



L2 



+ u ) dp 



(19) 
(20) 
(21) 

(22) 
(23) 



where 53 is the determinant of the 2+1-dimensional space-time giveii by (t, p, ip). Note that there is also another 
notion of energy in this space-time, namely that of the Tolman energy [2J, |25|. This defines the gravitationally active 
mass. 

In the Bogomolnyi-Prasad-Sommerfield (BPS) limit [13] given hy f3i = (32 = (3 — 2, 13^ — and with the choice 
q = g = 1 we have that = — such that it follows from (fTB]) that N{p) = 1. The remaining BPS equations are 



h' 



P' 



Ph 

T 

- 1 



/' = 

R!_ 
L 



Rf 



for the matter fields and 



£1. 
T 



-7 



2/l2p2 



2i?2/ 



2 -f 2 



-f {h^ - lY + if - lY 



L2 £2 

for the non-trivial metric function. The solutions fulfill an energy bound such that 



(n.m) 



27r(n + m) 



(24) 
(25) 

(26) 
(27) 



Note that in this limit the widths of the scalar cores pH,i become equal to the widths of the respective gauge field 
cores pw,i, i = 1, 2. 
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The binding energy per unit length of a (p,q)-string e\^'^^ can be defined as 



(28) 



Finally the (p,q)-string possesses magnetic fields in z-direction Bi = Bie^ and B2 = i?2ez with 



P' R' 



(29) 



where Bi and B2 are given in units of ^ . The magnetic fluxes then read 



<i>M,i = 27m , $M.2 = Svrm 



(30) 



and are obviously quantized. Hence, changing the winding numbers n and m changes the magnetic fluxes along the 
(p,q)-string. 



B. The geodesic equation 

The Lagrangian Cg describing geodesic motion of a test particle in the static cylindrically symmetric space-time 
(El) reads 



dx^ da;" 
dr dr 



N 



dr 



dr 



dr 



(31) 



where e = 0, 1 for massless or massive test particles, respectively and r is an affine parameter that corresponds to the 
proper time for massive test particles moving on time-like geodesies. The space-time has three Killing vectors ^ 

and ^ which lead to the following constants of motion: the energy E, the angular momentum Lz along the string 



axis (z-axis) and the momentum pz 



dr dr dr 



(32) 



Using the rescaling the constants of motion must be rescaled according to £^ — > E/{eir}i), pz — >■ Pz/{eirii), Lz 
Lz/{eir]if. We then find from ^ 



^ dr / \dT J \ dr 

Using the constants of motion wc find from ([31 

2 r^2 



e^-pI (dp 



1 (dp 

2\dT 



E^ 1 

2 2 



E' 1 



iV2 



7V2 ^ i2 



dr 



L2 



(33) 



(34) 



The left hand side of p4p is always positive and E"^ — e is & constant of motion. Following [2^ we can then rewrite 
this equation as 



(35) 



where 



VeSip) 



E' 1- 



(36) 



is the effective potential and £ = {E'^ — e)/2. Note that with this definition the effective potential becomes explicitly 
energy-dependent . 
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In the following, we would like to find t{p), (p{p) and z{p). For this, we rewrite the geodesic equation in the form 

d^ = ± (37) 

dz = ± —yj , (38) 

w(fa*-xgF-) 

The solution for each component can then be calculated as a function of p by using numerical integration methods. 



III. NUMERICAL RESULTS 



We have solved the set of differential equations ^TU\\ - p7|) numerically using the ODE solver COLSYS that uses a 
Newton-Raphson adaptive grid method [29]. The relative error of the solutions is on the order of 10"^'^ - 10~^°. Each 
component of the geodesic equation can then be integrated numerically by using the integrating function quad, i.e. a 
recursive adaptive Simpson quadrature in MATLAB with an absolute error tolerance 10~*. However the numerical 
profiles of the metric functions N{p) and L{p) must first be interpolated. This was done using a piecewise cubic 
Hermite interpolating polynomial, i.e. with pchip in MATLAB. With this procedure it is possible to obtain a smooth 
curve for the effective potential. 

In the following we will distinguish between bound orbits, escape orbits and terminating orbits. Note that when 
we talk about bound, escape and terminating orbits we are referring to the motion in the x-y-plane. The particles 
can, of course, move along the full z-axis from —oo to +oo for pz ^ 0. 

Bound orbits are orbits on which test particles move from a minimal value of p, /Omin > to a maximal value of p, 
Pmax < oo and back again. These orbits have hence two turning points with (dp/dr)^ = 0. On escape orbits, on the 
other hand, particles come from p = oo, reach a minimal value of p, p,nin > and move back to p = oo, which means 
that escape orbits have only one turning point with (dp/dr)^ — 0. Looking at ([55]) it is obvious that turning points 
are located at those p at which E — Vi^s (p) — 0. Finally, terminating orbits are orbits that end at the string axis p = 0. 

For all our calculations we have chosen g = 1 and g — ^■ 



A. Generalities 

Solutions to the model ([1]) have been extensively studied previously. The Table U summarizes the particular cases. 



Solution 


/3i,/32 
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Studied in 


Abelian-Higgs string in flat space-time 


/?i 




/33 =0 


7 = 


[2J 


Abelian-Higgs string in curved space-time 


/?i 


= /?2/0 




7/0 


[24], m 


(p,q)-string in flat space-time 


/3i 


= /32 = 2 


/33 /O 


7 = 


m 


(p,q)-string in curved space-time 


/?i 


= /32 = 2 


/?3 7^0 


7/0 


[15] 



TABLE L Known string solutions of tlie model ([T]). 

It has been observed in [23] that there are two types of solutions if one couples the Abelian-Higgs model minimally 
to gravity: string solutions and Melvin solutions which exist for the same values of the parameters in the model. 
These differ by their asymptotic behaviour of the metric functions at infinity. 



1. String solutions 



The string solution behaves like 

iV(p^Oo)=Ci, L(p oo) = C2P + C3, C2 > , 



(40) 
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where Ci, C2 and C3 are constants depending on n, m, 7 and /3i, z = 1, 2, 3. For /?3 = it has been found [2J, [23 
that ci > 1 for /3i = /32 = ;3 < 2, ci < 1 for /3i = /32 = ;3 > 2 and Ci = 1 in the BPS hmit fSi ^ P2 = P = 
A solution with the asymptotics (|40|) describes a conical space-time with deficit angle 6 given by 

5 = 2^(1 - C2) (41) 

In linear order the deficit angle S is given by the product of the coupling 7 and the inertial energy per unit length 
^(n+m) ^j-j-]^ J ^ gyc]^ |;j^g constant C2 = 1 for 7 = (or — Q) and C2 decreases for either 7 

or increasing. If the coupling 7 or the energy per unit length is too large then C2 < and the deficit angle 

S > 2tt. In this case the solution would have a singularity at a finite, parameter-dependent value of p = on with 
L{p — po) = 0, while N{po) stays finite. These solutions are the so-called super-massive string solutions [30| (or 
inverted string solutions) . We will not consider these kind of solutions in this paper and will always assume the deficit 
angle to be smaller than 27r. 

The "force" exerted on a test particle corresponds to the right hand side of 

•p, r'^"-pn„,^«V'. (42, 



dr2 \ J 

Note that for string solutions the effective potential tends asymptotically to a constant with Vcs{p — > 00) — > 
^ ^'^ 2c^^''"'"^' ^^'^ hence there is no force exerted on test particles far from the string. While the force associated 
to the angular momentum Lz is always repulsive, the total force close to the string can either be attractive or 
repulsive. Since — pj. > and iV > this depends on the sign of N' (see more details below) . 
For p <^ 1 the string solutions behave like 

7V(p < 1) - 1 + 0(p2) , i(p<l)^p. (43) 

Hence there is an infinite potential barrier at p = for test particles with non- vanishing angular momentum L^, i.e. 
these test particles can never reach the string axis at p = since their is no force to counterbalance the repulsive 
centrifugal force. On the other hand, for Lz = the effective potential tends to a constant Vosip — > 0) — ?► P2/2. 
Hence particles with — e < pi can reach the string axis. Since E^ > p\ these terminating orbits are only possible 
for massive test particles with e = 1. 



Infinitely thin cosmic strings The infinitely thin limit corresponds to the case where both the width of the scalar 
core as well as that of the gauge field core tend to zero. The string is hence a 1-dimensional object that can e.g. 
be described by the Nambu-Goto action. In this case the metric function N{p) = 1 (or some other constant that 
can be absorbed into the definition of t) and L{p) = C2p for p > 0. In this case, the only component in the force 
(1421) exerted on a particle is the repulsive angular momentum contribution. Hence, bound orbits are not possible 
in this space-time. This can also easily be understood when noting that the space-time of an infinitely thin cosmic 
string is locally flat @ and geodesies are just straight lines. The fact that bound orbits are possible in a finite width 
cosmic string space-time is related to the fact that close to the string axis the conical space-time is smoothed on scales 
comparable to the width of the string. The existence of bound orbits in "pure" cosmic string space-times [3l| is hence 
a new feature when considering cosmic strings with finite width. 



2. Melvm solutions 



The Melvin solutions exist for the same parameter values as the string solutions, but have a different asymptotic 
behaviour: 

N{p ^ 00) ^ aip2/3 , L{p -^00)-^ a2P~^l^ , (44) 

where again a\ and 02 are parameter dependent positive constants. This space-time is not asymptotically flat and the 
proper length of a curve with t = const., p = const., z = const and = — > 27r is s = 27ra2p~^/'^. This tends to zero 
for p 00. For the Melvin space-time with the asymptotic behaviour (l44l) the effective potential tends to infinity 
asymptotically with Vce{p — >■ 00) oc p^/^ for Lz ^ 0. Hence there is an infinite potential barrier at infinity for test 
particles with non-vanishing angular momentum, i.e. these particles can never reach infinity. This is related to the 
fact that the total force ([^^ on a test particle is always attractive at large p in Melvin space-times. For Lz — 0, the 
effective potential tends to E'^/2 for p — >■ 00. Hence, the asymptotic value of the effective potential is always larger 
than (for massive test particles) or equal to (for massless test particles) £. Massive test particles moving on radial 
geodesies can thus not reach infinity, while massless test particles have a turning point at infinity. 
For p ^ 1 the Melvin solutions behave like the string solutions P5)) . 
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B. Geodesic motion in (p,q)-string space-times: string solutions 

We will mainly discuss the geodesic motion in space-times with the asymptotic behaviour (j40[) since we believe 
this to be the physically relevant case. However, since the Melvin solution is a solution to the Abelian-Higgs model 
coupled minimally to gravity, we will also comment on this below. 



1. The effective potential 

The case /Ss = 0, /3i = /32 = /3 has been discussed for n = m = 1 in [2^. It was found that bound orbits are only 
possible for /? < 2 and for massive particles. In fact, in order to have bound orbits we need (at least) two turning 
points of the motion, i.e two intersection points between 14^ and £. Note that for £ finite and larger than the minimal 
value of the effective potential we will always have one intersection point for ^ due to the infinite potential barrier 
at small p such that escape orbits always exist. However, bound orbits are only possible if in addition the effective 
potential has local minima and maxima with ^^^^ = 0. At these local extrema we should then have 

LI N'{p) L{pf ■ ^ ^ 

Since - pi > Ojjy(p) > 0, L{p) > 0, L'{p) > this equation has only solutions for N'{p) > 0. For /Jg = it 
has been observed [26| that the metric function N{p) is either monotonically decreasing (for /? > 2) or monotonically 
increasing (for /3 < 2), while N{p) = 1 in the BPS limit (3 = 2. Hence the sign of N'{p) doesn't change and in 
particular, bound orbits are only possible for /3 < 2. In this case the energy-momentum part of the force (|^^ becomes 
attractive for (3 < 2, i.e. if the width of the scalar core is larger than the width of the gauge field core and can balance 
the repulsive part associated to the angular momentum. On the other hand for /3 = 2 (/3 > 2) the width of the scalar 
core is equal (smaller) than the width of the gauge field core. We observe that this leads to a vanishing (repulsive) 
energy-momentum part in the force (j42p and only escape orbits are possible. 

This is different when /Ss > 0. We will first discuss the case n = m = 1. The behaviour of the metric function 
N{p) of a (l,l)-string for 7 = 0.3 and different choices of (3i, i = 1,2,3 is shown in Fig. [TJ In all cases the blue 

dotted-dashed line corresponds to /Sa = and for cases (a), (b) and (c) the green solid line corresponds to /Sa ~ /Sg™'^'^^ 
with pI^^'^^^ = \/l3il32/2 the maximally allowed value for a given choice of j3i and /32- For (a) /?i = 1, = 2, (b) 
Pi — P2 — '2- and (c) /?! = /32 = 1 the increase of /Sa leads to an increase of the asymptotic value of N{p) for all choices 
of the i = 1,2,3. Hence, the increased binding between the p- and the q-string pronounces the effect already 
observed in the /Ja = limit. Note that while bound orbits are not possible in the BPS limit /3i = /32 = 2 for /3a = 
bound orbits do exist for /Sa > and /3i = /32 = 2 (which, of course, no longer corresponds to a BPS limit). For (d) 
/?! = 1.5, /32 — 6, (e) /?! = 2,/32 = 4.5 and (f) /3i = 2.25, /32 = 4 the metric function N{p) can have a local minimum 
if /3a < 132)1 i-G. if the binding between the strings is not too large. This is new as compared to the /3a = 

limit. We find that /3^"^(1.5,6) w 0.7, /3^"^(2,4.5) w 0.66, /3^"^ (2.25, 4) w 0.51. 

Obviously, for (a) /3i < 2, /32 = 2, (b) /3i = /32 = 2 and (c) /3i < 2, /32 < 2 the metric function N{p) increases 
monotonically, while for the other cases N{p) can first decrease from N{p = 0) = 1, have a local minimum at p = pmin 
with iV(/Omin) < 1 and then increase again to N{p ^ 1) > 1. This has important consequences for the shape of the 
effective potential as discussed below and can be understood as follows: consider the (l,l)-string to be a superposition 
of a (l,0)-string and a (0,l)-string. Now, for /3a = 0, these two strings do not interact. In this case, we know that for 
/3i > 2, j = 1, 2 the metric function N{p) would monotonically decrease, while for (3j < 2, j ~ 1,2 the metric function 
N{p) would monotonically increase. Superposing a string with j3i > 2 and one with /3j < 2 leads than to a metric 
function N{p) that first decreases and than increase again. Note that the opposite is not possible since the scalar core 
of a string with (3i> 2 \s smaller than that of a string with 13 j < 2. 

This behaviour of the metric function N{p) leads to the observation that the effective potential can have a negative 
minimum which for = is located exactly at Pmin- In the case (3^ = 0, the effective potential can have a local 
minimum for /3 < 2, but this will always be positive valued since N'{p) > means N{p) > 1. This is shown in Figl^] 
for a particle with E = 0.01, = 0.03 and Pz = 0. For cases (a), (b) and (c) the effective potential is positive for 
all values of p for non- vanishing or pz, while it can become negative for the other cases. In fact, the potential 
becomes positive everywhere for /3a ~ /^a'^'^^ • This will have influence on the existence of bound orbits as discussed 
below. In particular if the potential has a negative valued minimum as is e.g. the case for /3i = 1.5, /32 = 6 and 
/3a = 0.7, particles with f < 0, i.e. E^ < I can move on bound orbits. 

We have also investigated how the metric function N{p) changes when changing the winding numbers n, m and 
hence the magnetic fluxes along the string. Our results are shown in Figl3]for 7 = 0.3. 
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FIG. 1. The metric function N{p) of a (l,l)-string is shown for 7 = 0.3 and different choices of fii, P2 and Pz- 



We observe that the increase in the total magnetic flux along the string increases the asymptotic value of the metric 
function N{p) if N'{p) > 0. The qualitative features do not change. If a minimum of the metric function exists for 
n = m = 1 it exists for all choices of n and m (see Fig |3(d)| and if A^'(p) > for n = m = 1 this will be the same for 
other choices of n and m. 

Note that the profiles of the metric functions L{p) for all cases are similar to those for /Ja = 0. The deviation of 
L'(p ^ 1) — C2 from one determines the deficit angle of the space-time and depends on the inertial mass per unit 
length. This is shown for 7 = 0.3 and different choices of n, m and Pi in FigU) For solutions with /?i = (^2 the slope of 
L{p) at infinity decreases with increasing sum n + jn = p + q. This is natural since an increase in the windings leads 
to an increase in the mass per unit length and hence to an increase of the deficit angle. Moreover, for a given sum 
p + q the solutions with p = q have the biggest slope of L{p) at infinity. This is related to the fact that the solutions 
with equal winding are bound strongest (see also the results in jlsl). 



2. Classification of solutions 

The geodesies can be classified according to the test particles energy angular momentum Lz and momentum in 
the direction of the string axis Pz- Intersection points of £ with the effective potential, i.e. points where £ — Vcs 
correspond to turning points of the motion. The maximum and minimum of the effective potential determine the 
largest, respectively smallest possible value of £ for bound orbits. The effective potential is determined by the metric 
functions N and L as well as the constants of motion. Choosing /J^, i = 1,2,3, 7 and n, m we find the numerical 
profiles of iV and L. For a given (and e = or e = 1) there is an such that the value of £ is equal to the 
maximal value of the effective potential Vcs{p) and one E^ such that £ is equal to the minimal value of Vcs{p). In 
the former case, the corresponding orbit is an unstable circular orbit, while in the latter the orbit is a stable circular 
orbit. 

Defining p := E^ and v := we can then plot the domain of existence of bound orbits in the p-v-plane. Our 
results for massive particles with pz = are given in Fig.s[5l [5] for n = m — I. Fig IS] corresponds to the case of a 
(p,q)-string space-time with monotonically increasing N{p) and Fig. [5] to the case of a (p,q)-string space-time which 
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= 1.00, = 2.00, y = 0.30 



= 2.00, P^ = 2.00, Y = 0.30 



P^ = 1.00, P2 = 1.00, y= 0.30 



P3 = 0.00 
-P3 = 0.35 
-P, = 0.70 




0.14 
0.12 



0.06 
0.04 
0.02 



= 0.00 
= 0.50 
= 0.95 




0.14 
0.12 
0.1 

O.OB 

etf 
0.06 

0.04 



P3 = 0.00 
-P3 = 0.25 
-P =0.49 




/3i = 2, /32 = 2 



(b) 




0.14 
0.12 



"0.06 
0.04 
0.02 


-0.02 




(c) 



0.14p 
0.12 



"0.06 
0.04 
0.02 


-0.02, 



/3i = 2, /32 = 1 

P^ = 2.25, p^ = 4.00, Y = 0.30 



= 0.00 
= 0.45 
= 0.51 
= 0.60 



(d) 



/3i = 1.5, 132=6 



(e) 



I3i = 2, 132 = 4.5 



(f) 



I3i = 2.25, 132 



FIG. 2. The effective potential Vcff(p) in the space-time of a (l,l)-string is shown for 7 = 0.3 and different choices of /3i, /32 
and Ps. Here E = 0.01, = 0.03 and = 0. 



has N'{p) = at some non- vanishing, finite value of p. The blue dashed and solid black line from {p,i,0) to {p,3,i'3) and 
(/i2,0) to (/i3,z^3) , respectively, represent the choice of {E, L^, Pz) for stable and unstable circular orbits, respectively, 
and bound orbits exist in the colored domain between the two bounding curves. (^3,1^3) corresponds to the largest 
possible values of p, and v for bound orbits. Ml denotes the domain in the /u-i^-plane in which E is smaller than the 
minimum of the effective potential and hence there are no solutions to the geodesic equation. M4 denotes the domain 
in which £ is larger than the maximum of the effective potential and only escape orbits are possible. In M2 and M3 
on the other hand bound orbits are possible. In M2 £ is smaller than the asymptotic value of the effective potential, 
but larger than the minimum of V^ff and only bound orbits are possible. In M3 £ is larger than the asymptotic value 
of the effective potential but smaller than the maximum of V^s- Hence, in M3 there are bound orbits, but escape 
orbits are also possible. 

For /?i = /32 = 2 and /33 > we find that ^1 = 1 for all values of /33, while /i2 as well as {fi^^v^,) increase with 
increasing ji^. While for /3i = = 2, /33 = no bound orbits exist at all [2^, bound orbits are possible for = 0.1 
and the domain of existence of bound orbits in the /i-i/-plane is extending for increasing jS^ (compare the plots for 

= 0.1 and Pz — 0.75). The existence of bound orbits in the limit where Mn^i = M^^^i, i = 1, 2 is new as compared 
to the Pz —Q case. 

This, however, is not the only difference as compared to the space-time of an Abelian-Higgs string. As stated above 
we find that it is possible to have negative valued minima of the effective potential in (p,q)-string space-times. This 
leads to the observation that massive test particles with p, <1 can now move on bound orbits. This is a new feature as 
compared to the Pz = Q case, where we had to require that /i > 1. This means that test particles with less energy can 
move on bound orbits in (p,q)-string space-times as compared to the Pz = ^ case, which corresponds to the space-time 
of two non-interacting Abelian-Higgs strings. This is clearly seen in Fig. [S]for 7 = 0.35, Pi = 10, P2 = 3.6 and two 
different values of Pz- While for AfH,i > My^^i, i — 1,2 and in the Pz = ^ limit no bound orbits exist [1^ they exist 
in a small domain of the /i-i/-plane for sufficiently large Pz- The extension of the domain in the /i-i/-plane for which 
bound orbits exist increases with increasing Pz, i.e. the values of pi, 112 and {pz, 1^3) increase. 

The change of the /i-i^-plot of a (p,q)-string with 7 = 0.2, Pi = 8, P2 = 0.5, Pz = 0.99 resulting from the change 
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---p = 2,q = 2 




/3l = /32 = 4, : 



1.9 



P = 2,q = l 

p = 3,q = 1 





(d) 



/3i = 1, ^2 = 4, ^3 = 0.99 



(e) /3i = 2, ft = 4.5, ft = 1.45 (f) ft = 0.5, ft = 2, ft = 0.45 



FIG. 3. The metric function N[p) of a (p,q)-string for 7 = 0.3, different choices of i = 1,2,3 and different choices of 
(p,q)= (n,m). 
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(b) 
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(c) 



ft = ft = 0.5, ft = 0.24 
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= 2,q = 
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= 1,q = l 
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(d) 



Pi = 1, /32 = 4, ft = 0.99 (e) ft = 2, ft = 4.5, ft = 1.45 (f) ft = 0.5, ft = 2, ft = 0.45 



FIG. 4. The metric function of a (p,q)-string for 7 = 0.3, different choices of ft, i = 1, 2, 3 and different choices of 

(p,q)= (n,m). 
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of the winding numbers (p,q)= (n, m) and hence the change of the magnetic fluxes along the (p,q)-string are shown 
in FiglT] Here we concentrate on the case of a string with /Oh,i < Pw,i (the p-string) interacting with a string that 
has ph,2 > Pw,2 (the q-string). Increasing the winding of the p-string while keeping the winding of the q-string fixed 
shifts fj,i and /i2 to lower values, while the difference fi2 ~ /ii slightly increases with increasing p. Hence bound orbits 
are possible in a slightly bigger domain of the /^-^-plane and in particular test particles need less energy to be able to 
move on bound orbits when increasing the winding of the p-string. On the other hand, increasing the winding of the 
q-string while keeping the winding of the p-string fixed increases the value of /i2, while fii is nearly constant. Again, 
the domain of existence of bound orbits becomes larger when increasing the winding of the q-string. 

For pz 7^ 0, the qualitative features are the same. We observe, however, that the whole domain of existence of 
bound orbits shifts to larger values of fi when increasing p^. This is obviously related to the fact the the effective 
potential V^s is energy- dependent (see 

In contrast to massive test particles, we find that massless particles can only move on escape orbits. This is very 
similar to what has been observed in the ^3 = limit [261] and agrees with the result found in [s^l which states that for 
a general cosmic string space-time with topology x S massless test particles must move on g eodesies that escape 
to infinity in both directions, i.e. closed geodesies are not possible. The assumption made in [32| is that E must have 
positive Gaussian curvature. To show that S has positive Gaussian curvature in our case, we rewrite the metric ([S]) 
for massless particles (ds^ = 0) moving in a plane parallel to the x-y-p\ane as follows 

1 

dt^ = j^dp^ + J^dip^ = g^jdx'dx' , i^l,2 (46) 

where cjij is the so-called optical metric [s^ of which the spatial projection of geodesies of massless particles, i.e. light 
rays are geodesies, jjij is the metric of the above mentioned 2-manifold S and has Gaussian curvature K given by 

L' L" 

K = —N'N - —N^ - {N' f + NN" . (47) 

For /?! = = 2 and — (the BPS limit) we know that = 1 and the Gaussian curvature is obviously positive, 
away from the BPS limit one has to use the numerical solution and compute the curvature. We find that for most 
values of /3i, (32, P?, and 7 the Gaussian curvature is indeed positive and our result is in agreement with that of [s^ . 
However, if /3i and /32 are sufficiently large and sufficiently small, we find that K can become negative close to the 
string axis. Though the theorem of [32j is not applicable here, we nevertheless find that bound orbits do not exist. 



3. Examples of orbits 

In Figl5]we show how a massive test particle with E = 0.995, Lz = 0.022 and pz = 0.011 moves around a (l,l)-string 
with 7 = 0.35, Pi = 10, (32 = 3.6 and different choices of (i^. Note that the orbit is not planar due to the fact that the 
test particle has momentum in z-direction (see ([38])). The red and blue circles indicate the width of the scalar cores 
and the gauge field cores, while the dotted circles denote the minimal and maximal radius of the orbit. For our choice 
of parameters the width of the gauge field cores is larger than that of the scalar cores. We observe that the larger 
the closer the test particle moves around the string. For = 2.1, the orbit extends out to roughly four times 
the radius of the gauge field core, while for P^ — 2.2 the maximal radius is only roughly twice that of the gauge field 
core. Stating it differently: for smaller P^ the test particle moves mainly in the exterior vacuum region of the string, 
while for larger P^ it moves mainly close to or inside the string core, where the matter fields are non-trivial. As stated 
above, bound orbits are only possible in a limited domain of the /x-z/-plane. Test particles with values of E and Lz 
outside of this domain will not be able to move on a bound orbit around the string and will escape to infinity. An 
example of such an escape orbit of a massive test particle is shown in FiglHl Since the radii of the scalar and gauge 
field cores are very small in comparison to the extension of the orbit, we denote the core by a blue dot. The blue 
dashed line indicates the minimal radius of the orbit. We observe that the test particle comes from infinity, encircles 
the string core once and then moves again away to infinity for P^ = 2.15 and P^ — 2.48. For P^ = 2.96 the particle 
gets simply defiected by the string without encircling it. In fact, the deflection is decreasing for increasing P^. This 
can be explained by the fact that the energy per unit length and hence the deficit angle decreases with increasing P^. 

We have also studied how the orbits change when changing the winding numbers (p,q)= {n,m) and hence the 
magnetic fluxes. This is shown in Fig. [10] for the bound orbit of a massive test particle with E — 0.9931, Lz = 0.01 
and Pz = 0.015 in the space-time of a (p,q)-string with 7 = 0.2, Pi — 8, P2 — 0.5, P3 = 0.99. While for (p,q)= (1, 1) and 
(p,q)= (1, 2) the test particle moves close to the core of the string, it can extend considerably into the vacuum region 
for (p,q)= (2, 1). Apparently, the change of the winding of the q-string which has ph,2 > Pw,2 mainly influences the 
perihelion shift of the orbit, which increases with increasing winding. On the other hand the increase of the winding 
of the p-string which has p-w.i > Ph,2 allows the test particle to move further away from the string core. 
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The change of an escape orbit of a massless test particle with the change of the windings is shown in FiglTTJ While 
for (p,q)= (1, 1) and (p,q)= (1,3) the test particle gets simply deflected by the string it encircles the string before 
escaping to infinity for (p,q)= (3, 1). Apparently, the change of the winding of the q-string which has pH,2 > Pw,2 
influences the deflection only slightly. On the other hand the increase of the winding of the p-string which has 
Pw,! > /3h,2 leads to an encirclement of the string. 



C. Geodesic motion in Melvin space-times 



In Fig 112] we show the effective potential for a Melvin space-time with /3i = /?2 = 0.4, (3-^ = 0.18 and 7 = 0.57 
and test particle parameters E — 10, Pz ~ 5, Lz — 1.1. In comparison, we also give the effective potential of the 
corresponding string space-time. Close to the z-axis the effective potential of the Melvin space-time is equivalent to 
that of a string space-time. In contrast to the string space-time the effective potential in a Melvin space-time will 
always have a local minimum (see also (|^5|)) and Ves{p — > 00) — 00 for Lz 7^ 0. Hence there will be only bound 
orbits, while escape orbits do not exist. As already mentioned this is related to the fact that the space-time is not 
asymptotically flat and particles can never reach infinity. This is interesting since in string space-times a massless 
test particle will always escape from the string, while a massive test particle can move on a bound orbit provided 
both energy E and angular momentum Lz are not too large. In Melvin space-times test particles will always move on 
bound orbits. In Fig. [T3]we show the bound orbit for a massless test particle with E = 10, Lz = 1.1, Pz = 5 moving 
in a Melvin space-time with f3i — (^2 — 0.4, (3^ — 0.18, 7 — 0.57. 



IV. OBSERVABLES 



Since we believe the string space-time to be the physically relevant case, we compute all observables in the space-time 
with asymptotic behaviour (PH)) . 



A. Perihelion shift 



The perihelion shift of a bound orbit of a massive test particle (e = 1) can be calculated by the following expression 

S^e=i - 2 , ^ ^^"'^ ^ TTTJ - 2^ > (48) 



Mpy ( 



Pi 1 
T9~ TyTTT -L 



N(p}^ L(p) 



where p-min and Pmax are the minimal and the maximal radius of the bound orbit, respectively. The dependence of 
the perihelion shift 5tpe=i of a bound planar orbit of a massive test particle with E — 1.01, Lz = 0.02, = on 



the binding parameter [3^ is shown in Fig 14(a) For this particular case, the perihelion shift is negative which means 
that the test particle moves from the minimal to the maximal radius and returns back to the minimal radius under 
an angle less than 27r. In asymptotically flat black hol e sp ace-times and even asymptotically flat space-times of black 
holes pierced by infinitely thin cosmic strings (see p^.l23j|) the perihelion shift is positive. In a Schwarzschild-(Anti)- 
de Sitter black hole space-time the positive (negative) cosmological constant gives a positive (negative) contribution 
to the perihelion shift. Since all observations point to a positive cosmological constant, we would expect that the 
perihelion shift is positive for astrophysically relevant black hole solutions. Note that in the space-time of an infinitely 
thin cosmic string alone no bound orbits exist (see discussion above) and hence it makes no sense to calculate the 
perihelion shift. On the other hand, the presence of an infinitely thin cosmic string in black hole space-times enhances 
the (positive) perihelion shift ^] . The fact that the perihelion shift can become negative in the case of finite width 
cosmic strings is hence related to the fact that the space-time is that of a smoothed cone close to the string axis. In 
fact, the absolute value of the perihelion shift increases with increasing /33, i.e. for increasing (3^ the change of the 
(p coordinate from the first to the second minimal radius decreases. This can be understood when considering the 
infiuence of on the effective potential. In fact, the potential becomes steeper when increasing [3^ and hence the 
difference between the minimal and the maximal radius for fixed values of E, Lz and Pz decreases. Moreover, we 
observe that the perihelion shift for a (p,q)-string with A/jj,! > Afvv,i, i — 1,2 has the largest negative value, while a 
(p,q)-string with Mn.i < M^^^i, i — 1,2 has the smallest negative value. 

We have also studied whether the perihelion shift is always negative and find that it becomes positive for cosmic 
string space-times with coupling constants chosen such that the deficit angle is close to 2tt and the values of E and Lz 
are large, i.e. close to the boundary of the /z-i^-domain in which bound orbits exist (see Fig.s [S][7]). In this case, the 
difference between pmin and pniax is quite large and the test particle shows mainly in the vacuum region outside the 
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string. We find e.g. for a p-q-string space-time with /3i = /32 = 4, = 1.02, 7 = 0.499963 and resulting deficit angle 
S/{27t) « 0.9943 that the perihelion shift of a bound orbit of a massive test particle with E ~ 1.0018, = 0.007 and 
Pz = is positive and has value 6ipe=i ~ 2.3695 rad. 



B. Light deflection 



Is it very important for gravitational lensing to understand how massless test particles move on escape orbits. 
The deflection of light (massless test particle, i.e. e = 0) by a (p,q) string can be calculated by the following 
equation 

5^..o = r -jj, - n (49) 



where Pmin is the minimal radius of the escape orbit of the massless test particle. 

The dependence of the deflection S(pe=o of a planar escape orbit of a massless test particle with E = 1.04, = 0.28, 



Pz = on is shown in Fig 14(b) The light deflection decreases when increasing jS^. This is not surprising since 
the energy per unit length of the (p,q)-string and with it the deficit angle decrease with increasing /Ss. 

Moreover, we observe that the light deflection for a (p,q)-string with MH,i > Mw,i, i = 1, 2 has the smallest value, 
while a (p,q)-string with A^n.i < Myj^i, i = 1,2 has the largest value. This is related to the fact that if the scalar 
(gauge) fleld cores dominate strings tend to attract (repel) each other, hence lowering (increasing) the total energy 
per unit length as compared to the BPS limit. This leads to a decrease (increase) of the deficit angle. 



V. CONCLUSIONS 



In this paper we have studied test particle motion in the space-time of a cosmic superstring that consists of p D- 
strings and q F-strings. We have studied the asymptotically conical string space-time as well as the Melvin space-time 
that has vanishing circumference of a circle at infinity. We observe that the binding between the strings has important 
effects on the motion of test particles in string space-times. In the (33 = limit which corresponds to the space-time 
of two non-interacting Abelian-Higgs strings and is qualitatively similar to that studied in [26] massive test particles 
can only move on bound orbits if the scalar core width of the string is larger than that of the gauge field core. For 
> massive test particles can now move on bound orbits if the scalar core width is smaller than the gauge field 
core width and need less energy than in the /Ja = limit to do so. The perihelion shift can become negative due to 
the smoothed conical nature of the space-time close to the string axis and the absolute value of the perihelion shift 
increases with increasing /Sa . The fact that the perihelion shift can become negative seems to be a characteristic of the 
space-time of a finite width cosmic string that - to our knowledge - has not be noticed in any other astrophysically 
relevant space-time yet. Massless particles can only move on escape orbits and the deflection by the string decreases 
with increasing binding between the p- and the q-string. In Melvin space-times, on the other hand, massless and 
massive particles cannot escape to infinity and must move on bound orbits. 

The defiection of light by cosmic strings should be detectable. Though the identification of cosmic strings due to 
their gravitational lensing effects has been discussed extensively {33] no such cosmic string lens has been detected to 
date. Moreover, there might also be other sources of gravitational lensing and when identifying cosmic strings through 
gravitational lensing it has to be made sure that no "standard" matter distributions are the source of the lensing. 
On the other hand, the negative perihelion shift seems to be generic to finite width cosmic string space-times. To 
state it differently: if a negative perihelion shift would be observed this would be a strong evidence for the existence 
of cosmic strings. The main characteristic of cosmic superstrings is that they can form bound states and our field 
theoretical solutions describe such bound states. The fact that bound states can form so effectively alters the set of 
possible orbits considerably in comparison to standard field theoretical cosmic string models. In particular, the mass 
ratios 7 and /3i have an important impact. E.g. the perihelion shift can become positive or negative depending on 
the choice of these parameters and its absolute value can vary considerably. 
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FIG. 5. n-v plot for a (l,l)-string space-time with a monotonically increasing N{p). Here 7 = 0.3, Pi = (32 — 2 and ft — 0.1 
(left) and ft — 0.75 (right), respectively. The blue dashed and solid black line represent the choice of {E, Lz, Pz) for stable 
and unstable circular orbits, respectively. Here pz = 0. 




FIG. 6. The fi-v-plot for a (l,l)-string space-time which has N'(p) = at some non-vanishing, finite value of p. Here 7 = 0.35, 
Pi = 10, P2 = 3.6 and two different values of ft = 2.15 (left) and ft = 2.38 (right), respectively. The blue dashed and solid 
black line represent the choice of {E, Lz, Pz) for stable and unstable circular orbits, respectively. Here pz — 0. 
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FIG. 7. The /i-i/-plot for a (p,q)-string with 7 = 0.2, /3i = 8, /32 = 0.5, /Ja = 0.99 and different choices of (p,q)= (n, 
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3-d bound orbit 




X 

(c) /33 = 2.2 



FIG. 8. The bound orbit of a massive test particle with E — 0.995, Lz = 0.022 and pz — 0.011 in the space-time of a (l,l)-string 
with 7 = 0.3, /3i = 10.00, /32 ~ 3.6 and different choices of /Ja. The dotted circles denote the minimal and the maximal radius 
of the bound orbit, while the blue and the red circles indicate the radius of the gauge and scalar field cores, respectively. 




FIG. 9. The escape orbit of a massive test particle with = 0.02, E = 0.06 and pz — 0.05 in the space-time of a (l,l)-string 
with 7 — 0.35, /?! — 10, /32 = 3.6 and different choices of /?3. The blue dot denotes the core of the string, while the dashed blue 
line denotes the circle with minimal radius of the orbit, i.e. closest approach of the particle to the string. 
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3-d bound orbit 




(c) p = 1, g = 2 

FIG. 10. The bound orbit of a massive test particle with E — 0.9931, Lz = 0.01 and — 0.015 in the space-time of a 
(p,q)-string with 7 — 0.2, /3i — 8, 132 — 0.5, /3a — 0.99 and different choices of (p,q)= (n, m). The dotted circles denote the 
minimal and the maximal radius of the bound orbit, while the blue and the red circles indicate the radius of the gauge and 
scalar field cores, respectively. 
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FIG. 11. The escape orbit of a massless test particle with = 0.2236, E = 0.08 and pz — 0.05 in the space-time of a 
(p,q)-string with 7 — 0.2, /3i — 8, l32 = 0.5, /^a = 0.99 and different choices of (p,q)= (n, m). The blue dot denotes the core of 
the string, while the dotted line denotes the circle with minimal radius of the orbit, i.e. closest approach of the particle to the 
string. 
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E = 10.000, = 5.000, = 1.100, = = P3 = 0.18, y= 0.57 




FIG. 12. The effective potential 14ff (p) is shown for a Melvin solution (dashed) and in comparison for a string solution (solid) 
for E — 10, Pz — 5, Lz — 1.1 and coupling constants /3i = /32 = 0.4, /?3 — 0.18 and 7 = 0.57. The horizontal solid line denotes 
the value of £. 
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FIG. 13. The bound orbit of a massless test particle with E = 10, Lz — 1.1, Pz = 5 moving in a Melvin space-time with 
Pi = 132 — 0.4, — 0.18, 7 — 0.57. The blue and the red circle indicates the radius of the gauge and scalar field cores, 
respectively. 
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FIG. 14. The dependence of the perihehon shift of a planar bound orbit of a massive test particle with E = 1.01, = 0.02 
on and (left) and the dependence of the deflection angle of a planar escape orbit of a massless test particle with E = 1.04, 
= 0.28 on /?3 (right). 



